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Abstract 

The first measurement of decay-time-dependent CP asymmetries in the decay 
> J/if) Ag and an updated measurement of the ratio of branching fractions 
B{Bg^ J/ipKg)/B{B^^ J/il;Kg) are presented. The results are obtained using data 
corresponding to an integrated luminosity of 3.0 fb“^ of proton-proton collisions 
recorded with the LHCb detector at centre-of-mass energies of 7 and 8 TeV. The 
results on the CP asymmetries are 

Aaf (Rs —^ AV’-Rs) = 0.49 ± q'II (stat) ± 0.06 (syst) , 

Cdir (Rs —^ UV’-Rs) = —0.28 ± 0.41 (stat) ± 0.08 (syst) , 

5'mix (Rs —^ AV’-Rs) = —0.08 ± 0.40 (stat) ± 0.08 (syst) . 

The ratio B{B^^ J/ipKg)/B{B^^ J/ipK^) is measured to be 

0.0431 ± 0.0017 (stat) ± 0.0012 (syst) ± 0.0025 {fs/fd) , 

where the last uncertainty is due to the knowledge of the B^ and B^ production 
fractions. 
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1 Introduction 


In decays of neutral B mesons (where B stands for a B^ or B^ meson) to a final state 
accessible to both B and -B, the interference between the direct decay and the decay via 
oscillation leads to decay-time-dependent CP violation. Measurements of time-dependent 
CP asymmetries provide valuable tests of the flavour sector of the Standard Model (SM) 
and offer opportunities to search for signs of non-SM physics. A measurement of this 
asymmetry in the B^ —)■ J/if) Kg decay mode allows for a determination of the effective CP 
phase 

c|>f{B^^J/^PK^g)^c|>, + Acj),, ( 1 ) 

where (pd is the relative phase of the B^~B^ mixing amplitude and the tree-level decay 
process, and Acpd is a shift induced by the so-called penguin topologies, which are illustrated 
in Fig. In the Standard Model, (pd is equal to 2/3 (^, where /3 = arg(— 
is one of the angles of the unitarity triangle in the Cabibbo-Kobayashi-Maskawa (CKM) 
quark mixing matrix |^. The latest average of the Belle and BaBar measurements 
reads siiKp'f = 0.665 ± 0.020 [^, while the recently updated analysis from LHCb reports 
sin (j)f = 0.729 ± 0.035(stat) ± 0.022(syst) g. 

Forthcoming data from the LHC and KEK e’''e“ super B factory will lead to an 
unprecedented precision on the phase 0^®. To translate this into an equally precise 
determination of the CKM phase /3, it is essential to take into account the doubly Cabibbo- 
suppressed contributions from the penguin topologies, which lead to a value for A(j)d 
that might be as large as 0(1°) [^|^. By relying on approximate flavour symmetries, 
information on A(j)d can be obtained from measurements of CP asymmetries in decays 
where the penguin topologies are enhanced. The B^ —)■ J/'ip Kg mode is the most promising 
candidate for this task [^[^|^. 

Assuming no CP violation in mixing [^, the time-dependent CP asymmetry in 

—>■ J/ip Kg takes the form 


^ ^ r(g(t)^ J/^ATg) - F(i?0(t)^ 

F(i?o(f)^ J/V^iCO) + F(i?o(f)^ J/^|JK|) ’ 

_ Fmix sin (Am^ t) - Cgir cos (Am^ t) 
cosh (AFs f/2) -|- Map sinh (AF^ f/2) ’ 

where F(i?°(f)—)■ J/'ipKg) represents the time-dependent decay rate of the B^ meson into 
the J/'ijjKg hnal state, and Am^ = mu — and AF^ = Fl — Fh are, respectively, the 
mass and decay width difference between the heavy and light eigenstates of the B^ meson 
system. The 5° —)■ J/'i/j Kg CP observables are defined through the parameter 



A(B0^ J/PjK^g) 


( 4 ) 


in terms of the complex phase 0* associated with the B^-B^ mixing process and the ratio 
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Figure 1: Decay topologies contributing to the >■ J/ipKg channel: (left) tree diagram and 
(right) penguin diagram. 


of time-independent transition amplitudes as 


_ 2TZe[Xj/^x°] 
1 + 


Cdir = 


_ 1 ~ 


1 + 


_ 2Im[Xj/^Kl\ 


where Cdir and represent direct and mixing-induced CP violation, respectively. In 
the Standard Model (j)f^ = 2 arg(—V)^V)^). A recent analysis predicts 


( 6 ) 


^Ar (5°^ = 0.957 ± 0.061 , 

Cdir J/^Kl) = 0.003 ± 0.021 , 

J/i^Kl) = 0.29 ±0.20 . 


Similar expression for Eqs. ([^ and ([^ are obtained for the )■ J/'ipK^ decay by replacing 
s ^ d. The observable Map is not applicable in the measurement of —)■ J/tlj because 
it is assumed that AT^ = 0 [^. 

This paper presents the first measurement of the time-dependent CP asymmetries in 
5°—>■ J/'ipK^ decays, as well as an updated measurement of the ratio of time-integrated 
branching fractions 13(5°—?• J/'ip Kg)/B{B^ ^ J/'ijjKg). This ratio was hrst measured by 
the CDF collaboration [^, while the previously most precise measurement was reported 
by LHCb in Ref. 10 . The analysis is performed with a data sample corresponding to an 
integrated luminosity of 3.0 fb~^ of proton-proton (pp) collisions, recorded by the LHCb 
experiment at centre-of-mass energies of 7 TeV and 8 TeV in 2011 and 2012, respectively. 

The analysis proceeds in two steps. The first step, described in detail in Sec. consists 
of a multivariate selection of R —)■ J/'ipKg candidates. In the second step a maximum 
likelihood ht is performed to the selected data. The ht model includes a prominent 
B^ —)■ J/'ipKg component, which is used to improve the modelling of the B^ —)■ J/'ipK^ 
signal. In addition, the measurement of CP asymmetries associated with —)■ J/'ipK^ 
decays offers a validation of the likelihood method’s implementation. However, the stringent 
event selection necessary to isolate the 5° —)■ J/'0 77° candidates limits the precision on 
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these two CP observables. Dedicated and more precise measurements of the 5°—)- J/'ijjKg 
CP observables are therefore the subject of a separate publication [^. 

For a time-dependent measurement of CP violation it is essential to determine the 
initial flavour of the B candidate, i.e. whether it contained a 6 or a 6 quark at production. 
The method to achieve this is called flavour tagging, and is discussed in Sec. |4j The tagging 
information is combined with a description of the B mass and decay time distributions 
when performing the maximum likelihood ht, which is described in Sec. The three CP 
observables describing the 5° —)■ J/'0 decays and two CP observables describing the 

is derived from the ratio R of htted B^ —)■ J/%Ij to B^ —)■ J/%Ij event yields as 


>■ J/'ijjKg decays are obtained directly from the ht. The ratio of branching fractions 
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(7) 


where /sei is a correction factor for differences in selection efficiency between B^ —)■ iFg 

is the ratio of 5° to 5° meson 
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and 5°—)- J/ipK^ decays, and fs/fd = 0.259 ± 0.015 
hadronisation fractions. The study of systematic effects on the ratio R and the CP 
observables is presented in Sec. The main results for the branching ratio measurement 
are reported in Sec. [^and those for the CP observables in Sec. 


2 Detector and simulation 


The LHCb detector 13 14 is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing 6 or c quarks. The 
detector includes a high-precision tracking system consisting of a silicon-strip vertex 
detector surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4 Tm, and three stations 
of silicon-strip detectors and straw drift tubes placed downstream of the magnet. The 
tracking system provides a measurement of momentum, p, of charged particles with a 
relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200GeV/c. The 
minimum distance of a track to a primary vertex, the impact parameter, is measured 
with a resolution of (15 -|- 29/pt) M-ni, where px is the component of the momentum 
transverse to the beam, in GeV/c. Different types of charged hadrons are distinguished 
using information from two ring-imaging Cherenkov detectors. Photons, electrons and 
hadrons are identihed by a calorimeter system consisting of scintillating-pad and preshower 
detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons are identihed 
by a system composed of alternating layers of iron and multiwire proportional chambers. 


In the simulation, pp collisions are generated using Pythia 15 with a specihc LHCb 
conhguration [^. Decays of hadronic particles are described by EvtGen 17 , in which 
hnal-state radiation is generated using Photos 18 . The interaction of the generated 


particles with the detector, and its response, are implemented using the Geant4 toolkit 19 
as described in Ref. (201. 
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3 Event selection 


Candidate B —)■ J/'ipKg decays are considered in the J/il) —)■ and —)■ 7r’''7r“ 

final states. The event selection is based on an initial selection, followed by a two-stage 
multivariate analysis consisting of artificial neural network (NN) classifiers 


21 


3.1 Initial selection 

The online event selection is performed by a trigger, which consists of a hardware level, 
based on information from the calorimeter and muon systems, followed by a software 
level, which applies a full event reconstruction. The hardware trigger selects at least 
one muon with a transverse momentum px > 1-48 (1.76) GeV/c or two muons with 
\/PT(/ii) 7 »T(fi 2 ) > 1.3 (1.6) GeV/c in the 7 (8) TeV pp collisions. The software trigger 
consists of two stages. In the first stage, events are required to have either two oppositely 
charged muons with combined mass above 2.7 GeV/c^, or at least one muon or one high-px 
charged particle (px > 1-8 GeV/c) with an impact parameter larger than 100 pm with 
respect to all pp interaction vertices (PVs). In the second stage of the software trigger 
the tracks of two or more of the final-state particles are required to form a vertex that is 
significantly displaced from the PVs, and only events containing J/ijj —)■ candidates 

are retained. 

In the offline selection, J/ijj candidates are selected by requiring two muon tracks to 
form a good quality vertex and have an invariant mass in the range [3030, 3150] MeV/c^. 
This interval corresponds to about eight times the mass resolution at the J/'ij] mass 

and covers part of the J/ijj radiative tail. 

Decays of —)■ 7r'''7r“ are reconstructed in two different categories: the first involving 
K! mesons that decay early enough for the daughter pions to be reconstructed in the 
vertex detector; and the second containing 77° that decay later such that track segments 
of the pions cannot be formed in the vertex detector. These categories are referred to as 
long and downstream, respectively. Long candidates have better mass, momentum and 
vertex resolution than those in the downstream category. 

The two pion tracks of the long (downstream) candidates are required to form a 
good quality vertex and their combined invariant mass must be within 35(64) MeV/c^ of the 
known 77° mass [^. To remove contamination from A—;- pvr” decays, the reconstructed 
mass of the long (downstream) 77° candidates under the assumption that one of its 
daughter tracks is a proton is required to be more than 6 (10) MeV/c^ away from the known 
A mass 22 . The 77° decay vertex is required to be located downstream of the J/'ip decay 
vertex, i.e. it is required to have a positive flight distance. This removes approximately 50% 
of mis-reconstructed 7?°—)- J/077*(892)° background. The remaining 7?°—)- J/'077*(892)° 
background is heavily suppressed by the first stage of the multivariate selection described 
below. 

Gandidate B mesons are selected from combinations of J/^ and 77° candidates with 
mass in the range [5180, 5520] MeV/c^ and a decay time larger than 0.2 ps. The 

reconstructed mass and decay time are obtained from a kinematic fit 23 that constrains 
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the masses of the and 7T~^7r pairs to the known J/ip and masses 22 , respectively, 


and constrains the B candidate to originate from the PV. A good quality ht is required and 
the uncertainty on the B mass estimated by the kinematic ht must not exceed 30 MeV/c^. 
In the case that the event has multiple PVs, a clear separation of the J/'ip decay vertex 
from any of the other PVs in the event is required, and all combinations of B candidates 
and PVs that pass the selection are considered. 


3.2 Multivariate selection 

The hrst stage of the multivariate selection focuses on removing the mis-reconstructed 
)■ J/'ijjK*[892)^ background that survives the requirement on the flight distance. 
It only affects the subsample of candidates for which the is reconstructed in the long 
category. The NN is trained on simulated )■ J/'ipK^ (signal) and )■ J/'ll)K*(892)^ 
(background) data and only uses information associated with the reconstructed pions and 
Kg candidate. This includes decay time, mass, momentum, impact parameter and particle- 
identihcation properties. The requirement on the NN classiher’s output is optimised to 
retain 99% of the original signal candidates in simulation, with a background rejection on 
simulated )• J/'ipK*(892)^ candidates of 99.55%. This results in an estimated number 
of 18 ± 2 )■ J/'ijjK*[892)^ candidates in the long Kg data sample surviving this stage 

of the selection. Their yield is further reduced by the second NN classiher, and these 
candidates are therefore treated as combinatorial background in the remainder of the 
analysis. 

The second stage of the multivariate selection aims at reducing the combinatorial 
background to isolate the small B^ —>■ J/'ll) Kg signal. In contrast to the hrst NN, it is 
trained entirely on data, using the B^ —>■ J/i/ Kg signal as a representative of the signal 
features of the B/. —>■ J/'ip Kg decay. Candidates for the training sample are those populating 
the mass ranges [5180, 5340] MeV/c^ and [5390, 5520] MeV/c^, avoiding the B^ signal region. 
The signal and background weights for the training of the second NN are determined using 
the sPlot technique and obtained by performing an unbinned maximum likelihood 
ht to the B mass distribution of the candidates meeting the selection criteria on the hrst 
NN classiher’s output. The ht function is dehned as the sum of a B^ signal component 
and a combinatorial background where the parametrisation of the individual components 
matches that of the likelihood method used for the full CP analysis and is described in 
more detail in Sec. |5l 

Due to diherences in the distributions of the input variables of the NN, as well as 
diherent signal-to-background ratios, the second stage of the multivariate selection is 
performed separately for the B candidate samples containing long and downstream Kg 
candidates. The NN classihers use information on the candidate’s kinematic properties, 
vertex and track quality, impact parameter, particle identihcation information from the 
RICH and muon detectors, as well as global event properties like track and PV multiplicities. 
The variables that are used in the second NN are chosen to avoid correlations with the 
reconstructed B mass. 

Final selection requirements on the second stage NN classifier outputs are chosen to 
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Figure 2: Mass distribution of B candidates at different stages of the event selection for the (left) 
long Kg and (right) downstream Kg sample. The data sample after initial selection (red, +), 
after the first neural net (green, x) and after the second neural net (black, •) are shown. Overlaid 
are projections of the fit described in Sec.[^ Shown components are J/xjjKg (dark blue, 

dashed), B^^ J/ipKg (red, dotted) and combinatorial background (turquoise, dash-dotted). 


optimise the sensitivity to the signal nsing Ns/\/Ns + N-q as hgure of merit, where 
Ns and N^ are respectively the expected number of signal and background events in a 
±30 MeV/c^ mass range around the peak. After applying the hnal requirement on the 
NN classifier output associated with the long (downstream) Kg sample, the multivariate 
selection rejects, relative to the initial selection, 99.2% of the background in both samples 
while keeping 72.9% (58.3%) of the B^ signal. The lower selection efficiency on the 
downstream Kg sample is due to the worse signal-to-background ratio after the initial 
selection, which requires a more stringent requirement on the NN classifier output. The 
resulting J/ijj Kg mass distributions are illustrated in Fig. 

After applying the full selection, the long (downstream) B candidate can still be 
associated with more than one PV in about 1.5% (0.6%) of the events; in this case, one 
PVs is chosen at random. Likewise, about 0.24% (0.15%) of the selected events have 
multiple candidates sharing one or more tracks; in this case, one candidates is chosen at 
random. 


4 Flavour tagging 

At the LHC, b quarks are predominantly produced in bb pairs. When one of the two quarks 
hadronises to form the B meson decay of interest (“the signal 5”), the other b quark 
hadronises and decays independently. By exploiting this production mechanism, the signal 
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i?’s initial flavour is identifled by means of two classes of flavour-tagging algorithms. The 
opposite side (OS) taggers determine the flavour of the non-signal 6-hadron 25 while the 


same side kaon (SSK) tagger exploits the fact that the additional s (s) quark produced in 
the fragmentation of a (5°) meson often forms a {K~) meson 26 . 

These algorithms provide tag decisions qos and gssK, which take the value -|-1 (—1) in 
case the signal candidate is tagged as a i? (B) meson, and predictions r^os and rjssK for 
the probability of the tag to be incorrect. The latter is obtained using neural networks, 
which in the case of the OS taggers are trained on B~^ —)■ J/t!) decays, while for the 
SSK tagger simulated 5°—)- events are used. 

The mistag probability predicted by the tagging algorithms is calibrated in data to 
determine the true mistag probability cu, by using control samples of several flavour-specific 
B mesons decays. This calibration is performed individually for the OS and SSK tagging 
algorithms; for the latter, different calibration parameters are used to describe the B^ and 
B^g mesons. For all events with both an OS and SSK tag decision, a combined tag decision 
and mistag probability is derived as described in Ref. 


25 


The figure of merit for the optimisation of a tagging algorithm is the effective tagging 
efficiency, £eff = ^tag(l — 2a;)^ where £tag is the fraction of candidates with an assigned tag 
decision. In the long sample for the B^ —)■ J/'ip mode, the OS and SSK taggers yield 
an Seff of (2.93 ± 0.06)% and (0.97 ± 0.12)%, respectively, while the sample with both an 
OS and SSK tagging response gives an Sefr of (1.02 ± 0.10)%. In the respective downstream 
Kg sample, the OS and SSK taggers yield an Ses of (2.74 ± 0.11)% and (1.45 ± 0.15)%, 
respectively, while the sample with both an OS and SSK tagging response gives an of 
(0.48 ± 0.04)%. The combined Seff of all three overlapping samples for the 7?° —)■ J/'ipKg 
mode is measured to be (3.80 ± 0.18)% and (4.03 ± 0.16)% in the long and downstream 
Kg sample, respectively. 

In the B^ —)■ J/ipK^ mode, the main contribution is provided by the OS taggers, 
where the combined 5efr is measured to be (2.60 ± 0.05)% and (2.63 ± 0.05)% in the 
long and downstream Kg sample, respectively. Although the SSK tagging algorithm is 
specifically designed for B^ mesons, a small, but non-vanishing effective tagging efficiency of 
(0.064 ±0.009)% and (0.098 ±0.013)% in the long and downstream Kg sample, respectively, 
is also found for B^ mesons if the tag decision is reversed. This effect originates from 
same-side protons mis-identified as kaons, and kaons from the decay of K*(892)° mesons 
produced in correlation with the B^. Both tagged particles have a charge opposite to 
those of kaons produced in correlation with the and thus require the SSK tag decision 
to be inverted. Additionally, mis-identified pions carrying the same charge as the kaons 
correlated with the B^ dilute the effect described above. The SSK tagging response for 
B^ candidates is studied on 7?°—)- J/'ipK*(892)^ candidates using both data and simulated 
events. 
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Table 1: List of the observables describing the and systems that are included as Gaussian 
constraints to the likelihood fit. 


Parameter 

Valne 


Parameter 

Valne 


Arrid 

0.510 ±0.003 ps-^ 


Amg 

17.757 ±0.021 ps-i 

6 


AF, 

O 

1 


AF, 

0.081 ± 0.006 ps-i 

6 



1.520 ± 0.004 ps 



1.509 ± 0.004 ps 

6 



5 Likelihood fit 

The —)■ J/'ip Kg CP observables are determined from an unbinned maximum likelihood 
ht. The data is htted with a probability density function (PDF) dehned as the sum 
of a signal component, a signal component and a combinatorial backgronnd. In 
total it depends on seven observables. The PDF describes the reconstrncted B mass 
{jrij/^Kl ^ [5180,5520] MeV/c^), the decay time {t G [0.2,15] ps), and tagging responses 
gos and gssK- Additionally, it also depends on the per-candidate decay time nncertainty 
estimate 5t and mistag estimates and ?7 ssk- The long and downstream Kg samples 
are modelled using separate PDFs but htted simultaneously. The parameters common to 
both PDFs are the two )■ J/'ipKg and three 5° —>■ J/'i/jKg CP observables, as well as 
the observables describing the B^ and B^ systems that are listed in Table 


5.1 Mass PDF 

The mass shapes of the B —)■ J/'ipKg modes in both data and simulation exhibit non- 
Gaussian tails on both sides of their signal peaks due to hnal-state radiation, the detector 
resolntion and its dependence on the momenta of the hnal-state particles. Each signal 
shape is parametrised by a Hypatia fnnction 27 , whose tail parameters are taken from 
simnlation. The 5° component is constrained to have the same shape as the B^ PDF, 
bnt shifted by the B^-B^ mass difference, which is a free variable in the ht. The mass 
distribution of the combinatorial background is described by an exponential function. 


5.2 Decay time PDF 

The decay time distributions of the two signal components, T(t, gos, ^ssKlhos, ^ssk), need 
to be corrected for experimental effects originating from the detector response and the 
event selection. This is done by convolving them with a resolntion model, TZ{t\5t)i and 
combining the result with an acceptance function, £{t), to give the experimentally observed 
decay-time distribntion 

7~(t, gos, 5 'ssKlhos, hssK) X 77(t — dt^ x 8{t) . (8) 
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The resolution model has an individual width for each candidate, described by the per- 
candidate decay-time uncertainty estimate 6t provided by the kinematic fit introduced in 
Sec. A finite resolution reduces the amplitude of the oscillating terms in the decay-time 
distribution by a factor V = exp {—S‘fA'm‘^/2) [^,29 , and thereby affects the precision of 
the time-dependent CP observables. This effect is larger for the rapid oscillations 

than for the B^-B^ oscillations. The 5* estimates are calibrated using a separate sample 
of prompt J/ifj decays, which are produced directly at the PV and combined with random 
Kg candidates. This sample is obtained through the same event selection as described 
in Sec. except for the requirement on the decay time of the B candidates. The decay 
time distribution of the prompt J/'ijj mesons is modelled by the sum of three Gaussian 
functions sharing a common mean. For the long (downstream) Kg sample, this resolution 
model leads to an average dilution factor of {V) = 0.73 ± 0.13 (0.72 ± 0.04). 

The decay time distribution of the two signal components is affected by acceptance 
effects due to the decay-time bias induced by the trigger selection, the initial selection 
requirements and, most importantly, the NN classifier outputs. The shapes of the B^ and 

acceptances are assumed to be equal and modelled using cubic b-splines |^. The 
acceptance function is obtained directly from the data. The B^ — )• J/'ifj Kg decay time 
distribution is described by a single exponential, assuming AF^ = 0. The lifetime of 
the i?°, tbo = 1.520 ± 0.004 ps [^, is constrained in the ht using a Gaussian function 
whose mean is fixed to the known lifetime and whose width accounts for the experimental 
uncertainty. This allows the acceptance parameters to be directly evaluated in the fit to 
the data. 

The background decay-time distributions are modelled using two exponential functions, 
describing empirically a short-lived and a long-lived component. 


5.3 Likelihood fit 

The results are obtained from a simultaneous fit of the long and downstream Kg samples, 
using both the OS and SSK tagging information. In addition to the five CP observables, the 
nuisance parameters describing the mass (9 parameters), acceptance (12), background decay 
time (6) and event yields (18) are floated in the fit. The observables Am^, r^o, Arug, r^o 
and AFs, parametrising the B^ and B^ systems, and the effective B production asymmetries 
Aprod(-B°) and Aprod(-Bs) of the long and downstream Kg samples are constrained using 
Gaussian functions. The production asymmetries are dehned in terms of the B production 
cross-section <j{B) as Aprod(7?) = {(^{B) — a{B))/{a{B) + a{B)). The statistical and 
systematic uncertainties on the constrained parameters are added in quadrature and 
treated together; the correlation p(Fs, AF^) = —0.271 between the decay width and 
decay width difference of the B^ meson is also included. The effective B production 
asymmetries, specific to the data sample used in this analysis, are obtained by reweighting 
the results binned in B transverse momentum and pseudorapidity given in Ref. 
obtained values are listed in Table [U 

The likelihood ht is cross-checked using two independent implementations, and is 
validated with large sets of pseudoexperiments to thoroughly test several aspects of the 


31 . The 
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Figure 3: Decay time distribution of B candidates in the (left) long Kg and (right) downstream 
Kg sample. The fit projection is shown as solid black line. Shown components are B^^ J/ijjKg 
(dark blue, dashed), B^ -> J/ipKg (red, dotted) and combinatorial background (turquoise, 
dash-dotted). 


Table 2: Effective B production asymmetries specific to the data sample used in this analysis. 


Sample Mode 


Value 


Long ihg 

B' 

Downstream Kg 

B 

Long Kg 

B 

Downstream Kg 

5' 


-0.0117 ± 0.0057 (stat) ± 0.0013 (syst) 
-0.0095 ± 0.0051 (stat) ± 0.0013 (syst) 
-0.041 ±0.032 (stat) ±0.003 (syst) 
—0.022 ± 0.024 (stat) ± 0.003 (syst) 


analysis. These also include the use of stand-alone event generators that produce samples 
independently of the ht implementations. In addition, the ht model is tested on simulated 
data, with signal only and with both signal and background components present. The 
results from the ht to the full data sample are compared to those from various subsamples, 
and to those obtained from a weighted ht to the —)■ J/'^ Kg candidates only. All tests 
agree with the expectations and no biases in the ht are found. 

5.4 Fit results 

The results of the —)■ J/tlj Kg CP asymmetries are 

Cdir J/^Kl) = -0.028 ± 0.034 (stat) , 

(5° ^ J/^Kl) = 0.719 ± 0.034 (stat) , 
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where the uncertainties are statistical only. They are compatible with the BaBar 32 , 
Belle 33 and latest LHCb results. The results of the 5°—)- CP asymmetries 


are 


^Ar (5°^ = 0.49 ± ^,11 (stat) , 

Cdir (-5°—^ = —0.28 ± 0.41 (stat) , 

5'mix (-8° —^ = —0.08 ± 0.40 (stat) , 


where the uncertainties are statistical only, and the observed event yields are summarised 
in Table 1^ The ht projections for the mass and decay time distributions are shown in 
Figs. and 1^ respectively. The statistical correlations between the —>■ J/i^Kg CP 

observables are p(^Ar, C'dir) = —0.07, p(^Ar, = —0.01 and p(C'dir, S'mix) = —0.06. In 
addition, there is a 0(10%) correlation between ^af and the average decay width F^ and 
decay width difference AF^, and a 0(10%) correlation between Finix and the production 
asymmetries. The conhdence intervals for the three B^ —>■ J/'ipK^ CP asymmetries are 
also calculated with the Feldman-Cousins method 
with the point estimates given above. 


34,35 , which gives consistent results 


Table 3: Fitted yields from the unbinned maximum likelihood fit. The uncertainties are statistical 
only. 


Yield 

Long iPg 

Downstream 

5°^ J/ijKl 

27 801 ± 168 

51 351 ±231 


307 ± 20 

601 ± 30 

Combinatorial background 

658 ± 37 

2 852 ± 74 


6 Systematic uncertainties 

A number of systematic uncertainties affecting the determination of the 5° —)■ J/ipK^ 
CP observables and the ratio of event yields R are considered. The main sources of 
systematic uncertainty are due to assumptions for modelling the different components 
of the multivariate PDF. These uncertainties are estimated using large sets of simulated 
pseudoexperiments, in which the shapes and parameters of the individual PDF components 
are varied. In the generation of the pseudoexperiments, the values of the parameters are 
hxed to the ones obtained in the ht to the data. For each individual pseudoexperiment, 
the htted values of the CP observables and event yields are compared between the nominal 
ht and an alternative ht in which some of the shapes or nuisance parameters are varied. 
The resulting differences between the ht values form a Gaussian-like distribution. The 
mean and width of this distribution are added in quadrature and assigned as a systematic 
uncertainty. 
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Table 4: Summary of systematic uncertainties. 


Source 


C'dir 

Smix 

Long 

R X 10^ 

Downstream 
R X 10^ 

Mass modelling 

0.045 

0.009 

0.009 

15.5 

17.2 

Decay-time resolution 

0.038 

0.066 

0.070 

0.6 

0.3 

Decay-time acceptance 

0.022 

0.004 

0.004 

0.6 

0.5 

Tagging calibration 

0.002 

0.021 

0.023 

0.1 

0.2 

Mass resolution 

0.010 

0.005 

0.006 

12.6 

8.0 

Mass-time correlation 

0.003 

0.037 

0.036 

0.2 

0.1 

Total 

0.064 

0.079 

0.083 

20.0 

19.0 


Following this strategy, the systematic uncertainty due to the chosen mass model is 
evaluated by varying the Hypatia tail parameters within their uncertainties, replacing the 
signal model with a double Crystal Ball function 36 , and replacing the background model 


with a second-order Chebychev polynomial. The latter variation has the largest impact on 
the CP observables and yield ratio, and is used to assign a systematic uncertainty. 

The systematic uncertainty associated with the decay time resolution is evaluated by 
varying the dilution of the resolution model, through changes of the resolution parameters, 
and by comparing the nominal model with one that includes a scale offset in the calibration 
functions for the per-candidate decay time uncertainty estimates. The largest impact on 
the CP observables and yield ratio originates from the limited knowledge on the decay 
time resolution of the long sample. This forms the dominant systematic uncertainty 
to the —)■ J/V' Kl CP observables. 

Systematic effects due to the modelling of the decay time acceptance mainly affect 
^Ar, and are evaluated by varying the empirical model for S{t). 

The systematic uncertainty associated with the tagging calibration is obtained by 
comparing the nominal calibration with the largest and smallest effective tagging efficiency 
that can be obtained through changes of the calibration parameters within their respective 
uncertainties. 

The mass resolution is assumed to be identical for the and signal modes, but 
it could depend on the mass of the reconstructed B candidate. This effect is studied 
by multiplying the width of the B^ mass PDF by different scale factors, obtained by 
comparing B^ and B^ signal shapes in simulation. These variations mainly affect the ratio 
of event yields. 

Finally, a correlation between the reconstructed B mass and decay time resolution is 
observed in simulated data. The impact of neglecting this correlation in the ht to data is 
also evaluated with the simulated experiments. 

The total systematic uncertainty and its sources are summarised in Table 
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7 Branching ratio measurement 

The measured ratio of branching fractions is calculated from the event yields using Eq. ([^. 
The selection efficiencies and their ratio /sei are evaluated using simulated data. As 
the simulated data are generated with different values for the lifetime r^o, decay width 
difference AT^ and acceptance parameters compared to those measured in the collision 
data, correction factors are applied. This leads to a ratio of total selection efficiencies of 
/sei = 0.972 ± 0.029 for the long sample and fsei = 0.987 ± 0.040 for the downstream 
Kg samples. 

Combining the results in Table with the systematic uncertainties in Table yields 

R (long) = 0.01104 ± 0.00072 (stat) ± 0.00020 (syst) , 

R (downstream) = 0.01170 ± 0.00059 (stat) ± 0.00019 (syst) 

for the long and downstream Kg samples, respectively. A weighted average of the com¬ 
binations R X /sei for the long and downstream Kg samples is performed, assuming that 
they are uncorrelated measurements. The measured ratio of branching fractions is then 
given by 

b\bo^ j/IkI) ^ ^ ^ ^ • 

where the third uncertainty is due to the uncertainty in fs/fd- 

Combining the ratio of branching fractions with the known B^ —)■ J/'ipK^ branch¬ 
ing fraction B{B^^ J/'i/j K^) = (8.97 ± 0.35) x 10“^ [^, which accounts for the differ¬ 
ence in production rates for the B^B~ and B^B^ pairs at the T(4S') resonance, i.e. 
V{B^B~)/V{B'^B^) = 1.058 ± 0.024 j^, the B^^ J/'ijjKg branching fraction is 

B{Bl^ J/i,Kl) = 

[1.93 ± 0.08 (stat) ± 0.05 (syst) ± 0.11 {fs/fd) ± 0.07 {B{B^ J/^K^)] x 10“® , 

where the last uncertainty comes from the B^ —)■ J/'/ K^ branching fraction. 

8 Conclusion 

This paper presents the first measurement of the time-dependent CP violation observables 
in the decay B^^ J/'ijjKg and an updated measurement of its time-integrated branching 
fraction. Both measurements are performed using a data set corresponding to an integrated 
luminosity of 3.0 fb“^ of pp collisions recorded by the LHCb detector at centre-of-mass 
energies of 7 and 8 TeV. 

The results on the CP observables are 

•Aaf = 0-49 ± oil (stat) ± 0.06 (syst) , 

Cdir (Bs —t J/'ijjKg) = —0.28 ± 0.41 (stat) ± 0.08 (syst) , 

5'mix {B^ —)• J/'ipKg) = —0.08 ± 0.40 (stat) ± 0.08 (syst) . 
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The large statistical uncertainties on these results do not allow for a conclusive comparison 
with the predictions in Eq. 0 nor do they provide constraints on the shift parameter 
affecting CP measurements in 5°—)- J/ipK^. 

The ratio of time-integrated branching fractions is measured to be 


This result is the single most precise measurement of this quantity, and supersedes the 
previous LHCb measurement 
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